Observational Limits on Patchy Reionization: Implications for S-modes 
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The recent detection of secondary CMB anisotropy by the South Pole Telescope places a con- 
servative bound on temperature fluctuations from the optical depth-modulated Doppler effect of 
Taooo < -\/T3/iK at multipoles £ ~ 3000. This bound is the first empirical constraint on reionization 
optical depth fluctuations at arcminute scales, rsooo = O.OOf TaoooZ/iK, implying that these fluctu- 
ations are no more than a few percent of the mean. Optical depth modulation of the quadrupole 
source to polarization generates B-modes that are correspondingly bounded as -B3000 = 0.003 Tsooo. 
The maximal extrapolation to the I ~ 100 gravitational wave regime yields Bioo = O.IT3000 and 
remains in excess of gravitational lensing if the effective comoving size of the ionizing regions is 
i? ^ 80 Mpc. If patchy reionization is responsible for much of the observed arcminute scale temper- 
ature fluctuations, current bounds on Bioo already require R < 200 Mpc and can be expected to 
improve rapidly. Frequency separation of thermal Sunyaev-Zel'dovich contributions to the measured 
secondary anisotropy would also substantially improve the limits on optical depth fluctuations and 
_B-modes from reionization. 
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I. INTRODUCTION 

Recent observations by the South Pole Telescope 
(SPT) Q, and the Atacama Cosmology Telescope 
(ACT) [3| are ushering in a new era in which our un- 
derstanding of secondary anisotropy in the cosmic mi- 
crowave background (CMB) will be revolutionized. Sec- 
ondary anisotropy is generated after recombination by 
gravitational and scattering processes. It is thus more 
dependent on astrophysical processes than the primary 
anisotropy that has been so useful in determining fun- 
damental cosmological parameters. However, certain re- 
lationships between the various CMB secondary observ- 
ables can be used to scale out and, in principle, determine 
the unknown astrophysics. 

In this Brief Report, we discuss the example of these 
scaling relations provided by patchy reionization. Taken 
as an upper bound to account for contributions from 
other secondary effects, SPT measurements limit opti- 
cal depth fluctuations on arcminute scales. These same 
fluctuations generate B-mode polarization [J and so the 
implied limits at arcminute scales are relatively free of 
both cosmological and ionization model assumptions. 

By making a maximal extrapolation to the degree 
scales relevant for gravitational wave detection, we place 
upper limits on the contamination by patchy reioniza- 
tion i?-modes. Conversely, observational limits on degree 
scale B-modes constrain the ionization model, in particu- 
lar the size distribution of the ionized regions, when com- 
bined with arcminute scale temperature measurements. 



II. SCALING RELATIONS 

Thomson scattering of CMB photons off free electrons 
in linear velocity flows generates temperature fluctua- 
tions via the Doppler effect. The first order effect from 



the mean optical depth during reionization is highly sup- 
pressed on subhorizon scales and the dominant contri- 
butions on arcminute scales reflect optical depth modu- 
lation. Optical depth modulations can arise from linear 
density fluctuations (Ostriker-Vishniac effect), non-linear 
objects (kinetic Sunyaev-Zel'dovich effect), or ionization 
fluctuations (patchy reionization). 

Given a measurement or an upper limit on the tem- 
perature anisotropy due to this effect, one can constrain 
optical depth fluctuations and the corresponding effect 
on B-mode polarization. For notational convenience, we 
call the temperature power assigned to the modulated 
Doppler effect at a multipole of £ = 3000 
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Likewise, we use the general shorthand notation 
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for optical depth {X = r) and B-mode polarization {X = 
B) fluctuations; for i?-modes, s denotes the contribution 
from a particular source field. 

A conservative interpretation of the SPT detection of 
secondary anisotropy is that it places an upper limit of 
"^3^000 < 13 /xK2 at 95% CL. This limit assumes that all 
of the measured anisotropy is assigned to the modulated 
Doppler effect and considers the detection as an upper 
limit For reasonable cosmological models, one would 
expect the thermal Sunyaev-Zel'dovich fluctuations from 
unresolved clusters and groups to dominate the signal. 
Thus T^^ll < 5 AiK^ might serve as a more typical, al- 
beit model dependent, limit 0. For this reason, we will 
preserve the dependence of our main results on T^qqI ■ 

Optical depth fluctuations with a power spectrum CJ"^ 
modulate the Doppler effect from a velocity field with 



600 



400 



200 - 



T — I — I — r 




_i I I i_ 



_l I I L. 



_l I I L. 



10 

z 



15 



20 



At plausible redshifts for rcionization, z ^ 10, the pre- 
dicted quadrupole for flat ACDM niodels is a nearly con- 
stant Qrms = 18.2 ± 0.5 /uK, using the same data as for 
the 

'^rms predictions in Fig. ^ The uncertainty in Qrms 
is negligible compared with Wrms. 

Modulation destroys the symmetry that produces only 
_B-modcs from scalar perturbations, generating equal 
power in E and B polarization. When scaled to the tem- 
perature spectrum, the B-mode power spectrum from the 
modulated quadrupole becomes Q 
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FIG. 1: Predicted range of v-^^^i^z) (shading: 68% CL region; 
curves: 95% CL region) for flat ACDM (light gray) and non- 
flat quintessence (dark blue) models constrained by current 
data. 



rms I'l-ms to produce a temperature power spectrum [3| 
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below the coherence scale of the flows. The factor of 3 
here accounts for the line-of-sight nature of the Doppler 
effect. Note that for simplicity we assume optically 
thin conditions throughout and ignore the 10-20% effects 
from the finite mean opacity during rcionization. Given 
the linear theory matter power spectrum today, Pnn(fc): 
the rms velocity is computed as 
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where Di is the linear growth function for density fluc- 
tuations. Figure [T] shows the range of Wrms(-z) allowed by 
current constraints on flat ACDM models and non-flat 
quintessence models with arbitrary variations in the dark 
energy equation of state at z < 1.7, including measure- 
ments of the CMB, supernovae, baryon acoustic oscilla- 
tions, and the Hubble constant as described in Ref. 
Note that predictions for Umis at z > 2 have little uncer- 
tainty from cosmological parameters including variations 
in curvature and dark energy, with the possible excep- 
tion of the presence of a substantial component of the 
energy density in dark energy or massive neutrinos at 
high redshift. 

The same optical depth fluctuations modulate the gen- 
eration of polarization from the primordial quadrupole 
with rms Qrms. In the Sachs- Wolfe approximation [J|, 
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where the initial curvature spectrum normalized at the 
scale fco is A^(fc) = Asik/ko)"" , t] is conformal time with 
77* evaluated at recombination, and 

r f , r r[(3-n,)/2]F[(3 + n,)/2] 
^^^^"^^^'^r[(4-n.)/2]F[(9-..)/2]- 



Note that this scaling relation remains true regardless of 
whether the optical depth fluctuations are due to density 
effects like the Ostriker-Vishniac anisotropy @ or ioniza- 
tion effects from inhomogeneous rcionization. The only 
difference between these effects is the effective redshift at 
which Qrms/i'rms IS evaluated, with fj-ms dominating the 
variations. To maintain generality, we keep both Qims 
and Wrms in the relations but scale their values to z 10. 

Modulation of the e^"^ screening of the primary E- 
modes also generates i3-modes below the coherence scale 
of the i?-modes, i.e. the £ ^ 10'^ damping scale of the 
primary anisotropy. On these scales, the screening B- 
modes take the form 0] 
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where 
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-E-rms = 6.4 /iK for the maximum likelihood ACDM model 
and varies little with cosmological parameters. 

From Eq. ([3]), optical depth fluctuations are related to 
the temperature anisotropy TjlggQ as 
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Hence for the strict upper limit of Tgopg < ^/IZ^lK, the 
rms fluctuation in t at a few arcminutes is T3000 < 0.003, 
i.e. no more than a few percent of the mean optical depth 
(e.g. f = 0.10 ± 0.02 for general rcionization histories 
at 6 < z < 30 constrained by 5-year WMAP data [1]). 
Models that predict a percent optical depth rms or ~ 10% 
of the mean as a typical fluctuation are observationally 
unviable (cf. H). 

For S-modes from the modulated quadrupole [Eq. ([7])] , 
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To place upper limits on the B-mode contribution we 
begin by assuming all of the temperature power is due to 
ionization fluctuations since the B-mode signal from den- 
sity fluctuations is expected to be well below that of grav- 
itational lensing at all scales [l3| • For a variety of simple 
analytic models of patchy reionization that assume com- 
pletely ionized, spherical bubbles with a lognormal distri- 
bution of bubble radii [111, |T2| , the modulated quadrupole 
contribution to -B-modes can be approximated as a single 
peak in B^ that rises as £^ at £ ^ ^poak and falls as 

at £ > £peak, where ^poak ~ 2 X 104(i?eff/Mpc)-l . 

The effective ionized bubble radius i?cff is determined 
by both the volume-averaged radius of bubbles Ry and 
the lognormal width of the bubble size distribution crini?: 
i?cff = Rv exp(2.5crf„fl). 

Let us start by taking a single bubble size, i.e. a delta 
function in the distribution with crini? — > 0. Then the 
shape of the patchy reionization signal, shown in Fig. [21 
is approximately described as a sharp peak^ 



FIG. 2: Extrapolation of the inferred SPT limit on mod- 
ulated quadrupole B-modes (red arrow ai £ = 3000 with 
^3000 < 13 mK^, z ^ 10) to limits at ^ ~ 100 (black dots) 
using the sharp peak model with a single bubble size of 
iJeff = 200 Mpc (black dashed lines) or R^g = 50 Mpc (black 
dot-dashed lines). For RcS ^ 40 Mpc, we take a more conser- 
vative upper limit shown by the horizontal dotted line. For 
reference, solid blue curves show the contributions to iJf ex- 
pected from lensing and gravitational waves with tensor-to- 
scalar ratio r — 0.03 for the best fit fiat ACDM model. 
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The _B-mode rms at £ ~ 100 scales with i?off as 
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Combining quadrupole and E'-mode screening modula- 
tion yields 
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for 7 ^ i?off/Mpc < 200 and saturates outside this range, 
as shown in Fig. [31 

This assumption of a sharply peaked spectrum pro- 
vides a conservative extrapolation of upper limits from 
i - 3000 to £ ^ 100 if €pcak ^ 500 (or i?off Z 40 Mpc). 
In this case, realistic bubble size distributions would pro- 
duce a flatter spectrum and thus less power at ^ = 100 
(see 0, Fig. 7). 

In the opposite limit of i?cff ^ 40 Mpc, a flatter dis- 
tribution would produce more power at ^ ^ 100 than 
as the total _B-modes from reionization at £ = 3000. Thus a sharp peak and so we conservatively assume b[qq > 
upper limits on T^qqq place stringent bounds on patchy 
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reionization B-modes given predictions for Wims, Qrms, 
and -Erms- For the strict upper limit t!^qIq < vTS/iK, 
si?Jo<0.0lMK. 



III. LARGE ANGLE B-MODE LIMITS 

To relate the B-mode contributions at ^ = 3000 to the 
large angle regime at £ ^ 100 relevant for gravitational 
wave studies we require a model for the optical depth 
fluctuations. Note that the screening contributions die 
off as white noise for £ < £a ^ 300 given the acoustic 
scale and cannot generate large contributions at £ ^ 100 
01 . We thus need to relate the modulated quadrupole 
contributions between the two scales. 
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Using this model to extrapolate from the SPT 
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limit at £ ^ 3000 to £ ^ 100 produces an upper bound 
shown in Fig. [3las the shaded excluded region. 

Note that as the bubble size increases, the upper limit 
at £ = 100 saturates. This saturation is essentially model 
independent and corresponds to a maximal slope of £~'^ 
connecting the power at £ = 100 and 3000. This strict 
upper bound is given by 
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^ Ringing in £ produced by a delta function distribution is 
smoothed out by any more realistic bubble distribution and by 
the width in rcdshift of the reionization transition. 
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FIG. 3: Upper limits on the B-mode power from patchy reion- 
ization at ^ = 100 as a function of the effective radius of 
ionized bubbles, using T^qqq < 13 ^iK^ and assuming z « 10. 
The thin solid line extrapolates from the upper limit on 
to ^ = 100 using the scahng of Eq. (|14[l . and the thick solid 
line shows the extrapolation assuming a flat spectrum. The 
shaded region above both of these curves is excluded. The 
B-mode power at I = 100 from lensing in the best fit flat 
ACDM model is shown as a dashed hue. 



With ringp < y/TS fiK, this limit still allows more power 



than the lensing B-modes, -Biqq' ~ (2.68 ± 0.25) x 
10"'^ fiK'^ (with fractional errors approximately scaled as 
25a{nch^) [13). 

Conversely, observational limits on Bioo can be com- 

(v) 

bined with measurements of T^qqq to constrain the bub- 
ble size distribution. If the distribution is dominated by a 
single effective comoving radius, Eq. (jl4p provides an es- 
timate of the required precision to place an upper bound 
on this radius of 



•cff 



200(Bl?oV0.1Tiot)^/'Mpc 



Current bounds of S^^o < 0.1 ^lK^ (95% CL) from BI- 
CEP [ll] in fact exclude i?cff ^ 200 Mpc for the maximal 
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IV. DISCUSSION 



We have shown that the SPT detection of secondary 
temperature anisotropy interpreted as a limit on modu- 
lated Dopplcr contributions from scattering of < a/TS /iK 
provide the first empirical bounds on optical depth fluc- 
tuations during reionization at a scale of a few arcmin- 
utes. Models that produce fluctuations on these scales in 
excess of a few percent of the mean are no longer viable. 
This limit on optical depth variations in turn produces 
a nearly model- and cosmology-independent limit on the 
_B-modc polarization from reionization of < 0.01 /iK at 
similar scales. We expect both limits to improve dra- 
matically once more frequency channels of the data and 
larger regions of the sky have been analyzed so that the 
presumably dominant thermal Sunyaev-Zel'dovich con- 
tribution can be better separated. 

Such upper limits can be turned into constraints on 
degree scale _B-mode contamination of the gravitational 
wave signal. The maximal allowed B-modes at £ w 100 
can still exceed those from gravitational lensing, but only 
if the effective radius of the ionized regions is greater than 
^ 80 Mpc. In fact, current direct limits on degree scale 
B-mode polarization require ionized regions ^ 200 Mpc 
for the maximal allowed fluctuations at arcminute scales. 

Hence expected improvements in both small angle tem- 
perature measurements and large angle B-modc polariza- 
tion measurements should rapidly advance our empirical 
understanding of reionization. 
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